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STATISTICAL PROPERTIES OF FILTERED PSEUDORANDOM 
DIGITAL SEQUENCES FORMED FROM THE SUM 
OF MAXIMUM-LENGTH SEQUENCES 


INTRODUCTION 


In 1964 Shepertycki [l] suggested the use of filtered, pseudorandom 
digital sequences as the test signal in a root-mean-square (rms) system 
evaluation configuration. The digital sequence generator suggested by 
Shepertycki is the maximum -length type with linear feedback. 

Filtered maximum-length sequences have been used extensively for 
pseudorandom signal generation and similar applications [2] . The results 
were acceptable for low-pass filtered sequences when the ratio of the sequence 
digital-clock frequency to filter -cutoff frequency was approximately equal to 
the number of stages in the sequence-generating register. However, for 
higher ratio values, a definite skewing of the amplitude density function was 
observed by Gilson [3] . He reported that the skewing effect could be averted 
by using a digital sequence composed of the modulo-two sum of the outputs 
of two maximum-length digital-sequence generators. His paper consisted of 
the results of experiments with 9- and 11-stage generators. Approximately 
Gaussian distributed probability density functions were obtained for clock-to- 
filter cutoff ratios of up to 80 times the number of stages in the maximum- 
length sequence generator. 

Recent experiments have partially confirmed Gilson’s results; however, 
it was found experimentally and by computer simulation that sequence genera- 
tion by summing two maximum-length sequences does not always result in a 
nonskewed amplitude probability density function [4, pp. 42-44] . The computer 
simulation also revealed that although the amplitude distribution of the sum of 
two maximum-length sequences might be symmetrical, it is still not always 
approximately Gaussian; it is often platykurtic 1 as compared with the normal 
distribution. 


1. Peaked compared with normal distribution. 



Several pertinent analyses have been performed upon maximum -length 
sequences in the past. These include Gilson's experimental analysis, an 
analysis performed by Roberts and Davis [5], a computer simulation per- 
formed by White [6l , and an analysis by Lindhplm [7, pp. 569-5761 . . 

Gilson's paper indicated that the sum of the maximum -length sequences 
might provide analog noise with near Gaussian statistics when filtered but gave 
no analysis procedure. Roberts and Davis attempted to establish a mathemati- 
cal basis for predicting the statistics of filtered maximum-length sequences, 
but only mean and variance tests of the distribution were considered. White 
performed a computer simulation basically paralleling Gilson’s experimental 
efforts. White's results are inconclusive because of the small number of 
simulations attempted, but he indicated that a theoretical evaluation of the 
pseudorandom sequences would be of great value. Lindholm performed a 
theoretical evaluation of maximum -length sequences based upon the moments 
of weights of M-tuples from the sequence. He calculated the first three mo- 
ments of several maximum-length sequences. Wainberg and Wolf [8] have 
used Lindholm’ s approach in an analysis of M-tuples from maximum-length 
sequences and confirmed his results. 

The purpose of this study is to develop mathematical formulas to pro- 
vide the basis for the design of pseudorandom signals intended for applications 
requiring accurate knowledge of the statistics of the signals. 

The analysis approach involves generating the pseudorandom signals 
from a filtered hybrid-sum sequence rather than from a filtered single 
maximum-length sequence. The hybrid sequence is formed from the modulo- 
two sum of k maximum -length sequences and is an extension of the sum se- 
quences formed from two maximum-length sequences that Gilson evaluated. 

The basic reason for the analysis of hybrid-sum sequences is to establish a 
large group of sequences with good statistical properties. It will be shown 
that this can be accomplished much more efficiently using the hybrid-sum 
approach rather than forming the group strictly from maximum-length 
sequences. 

The need for information concerning the statistical properties of filtered 
pseudorandom digital sequences is critical because of the need to develop 
improved low-cost equipment for testing multichannel communications systems. 
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FORMATION OF HYBRID-SUM SEQUENCES 


Figure 1 is a block diagram of a hybrid-sum-sequence generator. The 
sequence is formed from the sum of k maximum-length linear digital sequence 
generators. 

The length of a hybrid-sum sequence is 
k n 

L = n ( 2 y -l) , (1) 

y=i 


where n 

y 


is the number of stages in the yth generator. 


The upper bound on the number of sequences available with k-sequence 
generators is 

, n 

k 9 y 

U = n . (2) 

n 

r=i y 


In equation (2) it is assumed that the maximum-length generators com- 
posing the hybrid-sum generator have fixed numbers of stages. If, however, 
there are a total of memory elements available which may be distributed 

in any ’’allowable” fashion among several maximum-length generators, the 
upper bound on the number of sequences available is 


U 


K 

l 


i=l 


k. 

i 

n 

y = l 


n . 



n . 
yi 


(3) 


where there are K ’’allowable” combinations of the N memory elements, and 

tfcl ^ 

k. is the number of separate generators for the i ’’allowable” combination. 

The condition, ’’allowable,” signifies that none of the sequence lengths of the 
maximum-length generators can have a common term. In other words, all 
the terms 
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N MAXIMUM -LENGTH GENERATORS 


w 

o 
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Figure 1. Hybrid-sum sequence generator. 



L 

y 


(4) 


- 2 Y -1 


for a particular hybrid-sum generator are relatively prime. Table 1 gives the 
factored values of L for values of n ranging from 3 to 34. 

y y 

n 

TABLE 1. FACTORIZATION OF L = 2 Y -1 INTO PRIMES 

y 


n 

2 Y -1 

2 3 

__ 

1 

= 

7 

2f 

__ 

1 

= 

3X5 

2 5 

- 

1 

= 

31 

2 6 


1 

■ = 

3x3x7 

2 7 


1 

= 

127 

2 8 

— 

1 

— 

3X5X17 , 

2 9 

_ 

1 

= 

7 x 73 

2 io 

- 

1 

= 

3X11X31 

2 11 



1 

'= ' 

23 x 89 

2 12 

_ 

1 

= 

3x3x5x7x13 

2 13 

— 

1 

= 

8191 

2 ^ 


1 

. = 

3 x 43 X 127 

2 15 

_ 

1 

= 

7 x 31 x 151 

2 16 

- 

1 

= 

3 x 5 x 17 x 257 

2 17 

— 

1 

= 

131071 

2®. 

- 

1 

= 

3 x 3 x 3 x 7 x 19 x 73 

2 19 

_ 

1 

= 

524287 

2 20 

- 

1 

= 

3 x 5 x 5 x ll x 31 x 41 

2 21 

— 

1 

= 

7 x 7 x 127 x 337 

2 22 

- 

1 

= 

3 x 23 x 89 x 683 

2 23 

_ 

1 

= 

47 x 178481 

2 24 

— 

1 

- 

3 x 3 x 5 x 7 x 13 x 17 x 241 

2 25 

- 

1 

= 

31 x 601 x 1801 

2 26 

- 

1 


3 x 2731 X 8191 




TABLE 1. (Concluded) 


n 

y 

n 

2 Y -1 

27 

2 27 - 

■= 7 x 73 x 262657 

28 

2 28 - 

L = 3x5x29x43x113x127 

29 

2 29 - 

L = 223 x 1103 x 2089 

30 

2 30 - 

L = 3 x 3 x 7 X 1 x 31 X 151 x 331 

31 

2 31 - 

= 2147483647 

32 

2 32 _ 

= 3 x 5 x 17 x 257 x 65537 

33 

2 33 - 

L = 7 x 23 x 89 x 599479 

34 

2 34 - 

L = 3 x 43691 x 131071 


Table 1 can be used to select the number of stages n^ for each of the component 
maximum-length sequence generators of a hybrid-sum sequence generator. 

MODEL OF FILTERING PROCESS 

The analysis of the statistical properties of low-pass filtered hybrid- 
sum sequences is based upon a filter with an impulse response of the "brick- 
wall" type as shown in Figure 2. 



0 m b 


Figure 2. Impulse response of assumed filter. 
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The brick-wall approximation to a filter’s impulse response is equiva- 
lent to representing the impulse response as the sum of unit-step functions 

h(t) = u(t) - u(t - MT b ) , (5) 

where u(t) is the unit-step function. 

The transfer function of the assumed filter is 


( gdMTjA -ju>MT b 

— « — / ^ 

... „ . . 1 e • (6) 

-co — 

2 

This filter impulse response shown in Figure 2 serves as an approxi- 
mation to the impulse response of filters commonly used in the baseband of 
communications systems. The primary justification for using the approxima- 
tion in the analysis is the. degree of correlation between the results of the 
analysis and experimental data obtained by filtering sequences [4, pp. 45-125]. 

Using the brick-wall impulse response approximation, the filtering 
process is equivalent to forming a signal by adding digits of the sequence when 
the digits are ± 1. If the central-limit theorem is applied, the output signal 
density function tends toward a Gaussian density function if the sequence is 
random and the number of digits added is large. This is true because the 
individual digits forming the sum are independent for a random sequence [9] . 

The digits of a hybrid-sum sequence are not actually independent, but 
in many cases they approximate the random sequence well enough to approach 
a Gaussian distribution when filtered. 


STATISTICAL PROPERTIES 

An analysis of some of the statistical properties of hybrid-sum 
sequences is now presented. If x(t) is a digital sequence, its autocorrelation 
function is 
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( 7 ) 


R (t) 

XX 


lim 1^ 
T— co T 


T/2 

/ x(t)x(t+T)dt 
-T/2 


If x(t) is random, the autocorrelation function will have the form shown in 
Figure 3. 



Figure 3. Autocorrelation of a random sequence, x(t) . 

If x(t) is a maximum -length sequence (sequence length equals L) an appropri- 
ate autocprrelation function is 


R M = \ f x(t)x(t + T)dt , (8) 

T -T/2 


where T equals LT, , and T, is a bit period, 
b b 

Now assume that sequence x(t) is generated by summing (modulo-two) 
k maximum-length sequences, 

k 

x(t) = Y) x .(t) (modulo-two) (9) 

i=l 1 
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when x(t) is from { 1, 0} , or 


x(t) = n x.(t) 
i=l 1 


( 10 ) 


when x(t) is from {-1, +l} . 

The autocorrelation function from equation (8) becomes 


R (t) 

xx 


T/2 k k 

f n X (t) n X ( t + r ) dt 
T ■ y _ y 


T J , y 

-T/2. -7=1 7=1 


( 11 ) 


or 


R (t) 

xx 


T/2 k 

T / n X (t)x (t + T)dt 
-T/2 y=l r y 


( 12 ) 


When binary digits from { -1, +l} are used, the equivalent of the shift-and-add 
property becomes the shift-and^multiply property, and if 


t + qLy for 7 = 1, 2, 3, 


k, q is an integer, and r + 0, 


then 


= X ^ ft + T ^( T ^ 
yy y y 

Substituting equation ( 13) into equation ( 12) , 


(13) 


T/2 k 

R (t) = — f II x (t + r )dt 

XX T J . *v v -v 7 


-T/2 7=1 


y 


(14) 
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where each is some function of T. Since the sequence lengths have no 

common terms, all possible phase arrangements of the constituent sequences 
are cycled during the total sequence length T. This characteristic of the 
product sequence can be used to change the form of equation ( 14) to 


R (T*qL ) 
xx 7 


O 

-{III 

T 1 =1 V 1 V 1 


k k 

Z n X (t) 

T, =1 *y=l y 
k 


i X, x i ( V X: X 2 ( V ••• X. \ {t k> < 15 > 


V 1 


V 1 


V 1 


where 


r = n t 


r=i 


y 


(16) 


In each sequence there is now one more -1 than +1 and 


R ( T *0) = 

xx L 


(17) 


If r equals 0, there is no shift and 


R (t=0) ■= 1 

XX 


(18) 


For the case where r equals qL^ for 1< j< k, and q is an integer, then 


R (T=qL.) 
xx y 


T/2 k 

T / n *(t)x„(t + T). 


-T/2 7=1 

T^j 


y y 


(19) 
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From the shift and multiply property, 


R (T=qL.) 
xx ' J 


1_ 

T 


T/2 k 

f n x (t+r )dt 
-T/2 ..y=l 7 7 

y*) 


or 


L L 
1 2 


k k 


■ l Z Z ••• Z n 

T =1 T =1 T,=l 7=1 7 

12 K 

7^1 


( 20 ) 


= t i 


T=1 
] 


k 

n 

7=1 

7^j 


XJ 

i 

T =1 
7 


x (T ) 
7 7 


R xx (T=qL j ) " T ( - 1)k ' 1 


( 21 ) 


Figure 4 gives an example of the autocorrelation function as a function of r for 
even and odd values of k. The scale is exaggerated for illustration purposes, 
and the figure demonstrates that the autocorrelation of a repetitive function 
is repetitive. 

An analysis of the randomness properties of the filtered sequence can 
be based upon the calculation of the moments of the resulting waveform. 

The filter time output, r(t), is 


00 

r(t) = f x(r)h(t-r)dr . (22.) 

-oo . 
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With the assumed filter impulse response, 


M 

r(t) « X x ( t_1T b ) 
i-1 

where x(t) is the digital sequence. Equation (22) states that r(t) is equal to 
the sum of the last M digits that occurred before t, and the study of the statis- 
tical properties of the filter output is based on the weights of the sequence 
M-tuples. ' 



Figure 4. Examples of the autocorrelation function of 
sum sequences from k generators . 

th 

The i M-tuple weight of the sum sequence is 



where 



k 

n 

7=1 


i+j»T 


( 24 ) 


th 

Equation (24) equates the state of the (i+j) bit of the sum sequence 

th 

to the product (equivalent to modulo-two addition) of the (i+j) member of 
the k component sequences, 

The statistics of the weight functions can be investigated by calculating 

th 

the moments of the weight distribution. The p moment is 

S p =l r> s P • (25) 


Substituting equations (22) and (23) 

into equation (25) yields 



L-l 

( M-l 

k 

\ 


o p _ I 
S " L 

l 

i =0 

\ 3=0 

n 

7=1 

x. + . P . 

1 + 3 , y/ 

(26) 

The first moment is 





L-l 

M-l 

k 



sl = i 

Z 

i =0 

1: 

j=o 

n j 

7=1 

c i + j.r 

(27) 


M-l 

L-l 

k 



1 

L 

l 

l 

i= 0 

n x. . 
y=l ^ 
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M-l k / y 

- T. Z " l 


L -1 


j=o y=i 


i=0 


i + j»r 


(27) 
(Con. ) 


! (-D k 


The first moment does not depend on the particular sequences comprising the 

sum sequence, and if .M«L, then is approximately 0. The first moment 
indicates the mean of the weight distribution. 

The second moment is 


„ . L-l / M-l k 

» 2 = i ill » 

i=0 \ j=0 7=1 


x. 


i+j,7 


, L-l / M-2 M-l k 

I l I M + 2 l l n x .. x. 


i=0 


j=0 £=j+l 7=1 


i+j,7 i+£»7 


(28) 


M-2 M-l L-l k 


M + T Z Z Z n x. . x 

L . i . , i+J.7 i+S,7 


j=0 £=j+l i=0 7«1 


Now 


X X — X 

i+j,7 i+£,7 i + z»r 


(29) 


if M < L for all y, and z is some function of j, £, and y from the shift-and- 
7 

multiply property of maximum-length sequences. 


14 



Equation (28) becomes 




M-2 

l 


i=0 


M-l 

■l 


E=j+1 



( 30 ) 


Now 



and 


S 2 =M+-^- t ( M-l)M ] (32) 

L 

if M<L for all y. The second-central moment (the second moment about the 

y 

mean-of-the-weight distribution) is 



or 





(M-l) 



where M<L for all y. 

V 


(33) 


(34) 
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From equation (34), the second-central moment does not depend on 
the particular sequences comprising the sum sequence but does depend 
upon the number of sequences summed. 

As an example, consider the case where k equals 1, n^ equals 5, and 
L equals 31: 



The maximum second-central moment is for M approximately equal to 15, 

2 

where S c is approximately 8. 

Asa second example, consider the sequence generated by modulo-two. 
summing of two maximum-length sequences. The two sequences are from 
generators with two and three stages: 

n = 2, n = 3, L = 21 
1 z 

and (36) 



if M<3. 

The second-central moment is a measure of central tendency. Observe 
from the two examples that for these particular sequence generators the 
second-central moment for the summed sequence rs larger than the second- 
central moment of the single maximum -length sequence for 1<M< for ally. 

The generators in the two examples each require five register stages (case 1: 
n^ equal to 5; case 2: n^ equal to 2, and equal to 3) , and the sequences are 

approximately the same length. 
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The third-central moment of the weight function gives an indication of 
the relative skewness of the graphical representation of the function probability 
density. The third moment is 


Now 


, L-l / M-l k 

s 3 = i ill n x . 

L i = 0 V j=o y=i 1+J ’ Y 


(37) 


M-l k 

Z n x 

j=o y=l 


i + j,y 


M-l k 

(3M-2) l II x 

j=0 -pi ‘ J ’ Y 


M-3 M-2 M-l k 

• 3! Y v y ii x x , 

£o >=i 1+V - Y 1+V ' Y 


X. 


and 


(38) 


_ M-l L-l k 

g3 (3M-2) £ £ „ x 

L tfo V=1 1+) ’ Y 


+ M 

L 


Using the relations 


M-3 M-2 M-l L-l k 

ll l l n x. . x. x 

A v=j+l £=v+l i=0 y=l 1+J ’ Y 1+V ’ Y 1+ ^ Y 


(39) 


L-l k 
Z n x 
i=0 y=l 


i+j.y 



i+j.y 


(40) 
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and 


k 

n 

y=l 


X X X 

i+j,y i+v,y i+£,y 


n x. ■. x.,, 

7=1 i+j.r i + 0,r 


(41) 


where 0 is a function of v, f, andy, and M< L for ally, equation (41) 
becomes y 


S 3 = ~ (3M-2) (M) (-l) k 

Xj 


r . Z X E Z 11 Vj y X i+0 y • (42) 

j= 0 V=j+1 £=v+l i=0 . 7=1 hy 9,y 


M-3 M-2 M-l L-l k 


If 


0 = j 


(43) 


L -1 

l 

i=0 


x.j . x. . = L 

i+j,7 i+0,y i 


(44) 


and if 


0*3 » 


(45) 


then 
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(46) 


L -1 


i =0 


x i+j,r x i+<t>,y - 1 


where M < L for all 7 . 
7 . 


There are ^^combinations of j, v, and £ in the latter part of equation (42). 

Let C be the number of these combinations for which the set o> (with w 

w 

members) of the k component sequences satisfies 


X i+j,7 X i+v ,7 X i+£,7 


(47) 


for all i. 


The third moment becomes 


S 3 = i (3M-2)(M)(-l) k + ^ l 

0 ) = 1 


w 

c n l (- 1 ) 

co y 

7=1 


k-w 


,(48) 


where M<L for all 7 . 

7 

If equation (47) is satisfied for a particular component-maximum-length 
sequence, then that sequence must satisfy the trinomial equation 


+ A ^ + .A ^ - V = 


(49) 


where A is the generator A matrix. The minimum polynomial that describes 
the sequence is . 
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(50) 


n n , 

A T + « t A 1 '- 1 •+ 


... + a 

n-1 


+ I 


0 


and is called the characteristic equation. If the A matrix satisfies both equa- 
tions (49) and (50) , the polynomial on the left of equation (50) is a factor of 
the polynomial on the left of equation (49) . . If the maximum-length sequences 
under consideration have characteristic equations that are factors of no 
common trinomials of power equal to or less than M-l, then equation (48) 
becomes 


S 3 = ~ (3M-2) (M) (-l) k + H 


0 - i *, 


(-D + I bl(j) 

y+1 Y Y 


k-1 


(51) 


if M < L for ally, where B is the number of trinomials of power less than_or 

y y th 

equal to M-l that contain the y -sequence characteristic polynomial as. a factor. 
If M«L, then 


+ 3 1 


n l 


y=l 


VV ( - l) 


k-1 


(52) 


7=1 


3 

where S is the third-central moment and M<L for all y. 
c y 

From equation (52) , trinomial factors tend to increase the third- 
central moment. To minimize the third-central moment, sequences should 
then be added that have characteristic equations that are factors of few trino- 
mials of power less than or equal to M-l. 

The fourth moment is 


/ 
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. . L-l'/M-l k 

s 4 = i III n x 

^ i=0 V j=0 7=1 J,T 


( 53 ) 


or 


. n L-l 

s 4 - i E 


i=0 


M-2 M-l k 

M + 12M(M-1) + 2M(3!) J\ T, n x. . x. 

j=0 v=j+l y=l J,K 


M-4 M-3 M-2 M-l k 

+ 41 E E E E n VvVvVtr x Wr 

j=0 v=j+l £=v+l ^£+1 7=1 J,T ,y e,y 1 1,7 


(54) 


Now using the shift-and-multiply property of maximum-length sequences, and 
assuming M<L for ally, 

y 


A /n ,_i o y M-2 M-l k / ' y_l 

S 4 = M + 12M(M-l) + — Yj Z 11 ( Z 


x. 


j=0 v=j+l y=l\ i=0 


i+x,y 


. , M-4 M-3 M-2 M-l k / y- 

+ ii e E E E n E x. 

L j=0 v=i+l £=v + l *.£+1 t=i\ 1=0 1+Z 1' Y 1+ V Y 


/-I 


(55) 


There are combinations Of j, v, £ , and f in the latter part of 


equation (55) . Let D be the number of these combinations for which the set 


co 


co (with w members) of the k component sequences satisfies 


X — x X X 

i+j,y i+v,y i+£,y i+f,y 


( 56 ) 


for all i . 
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Then 


s 4 = M + 12M(M-l) + -1^5- (M-l) t-l) k 


T l 


IX> = 1 


W 

d n l (- 1 ) 
w , y 
y=i 


k-w 


( 57 ) 


If equation (56) holds, the particular sequence satisfies the quadrinomial 
equation 


i-v i+t i-f 
x J + x J b + x J +1=0 , 


(58) 


where j-v ^ M-l. 


From the properties of the characteristic equation for maximum-length 
sequences, if equation (56) holds, the characteristic equation will be a factor 
of the quadrinomial in equation (58) . If none of the characteristic polynomials 
of the constituent maximum-length sequences are factors of common quadri- 
nomials of powers less than or equal to M-l, then equation (57) reduces to 


S 4 = M + 12M(M-1) + 12 _ M ■ (M-l) ( -l) k 


+ 11 
L 




E 


y 


(-D 1 


t E L (.1) 

y=l 77 


k-1 


(59) 
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where E is the number of quadrinomials of power less than or equal to M-l 


y 


th 


that have the 7 characteristic equation of a factor. 
If M«L, 


l E y L y (- 1 ' 

S 4 W M + 12M (M-l) + 4! 


k-1 


n l 

y=l ^ 


4 

where S is the fourth central moment, 
c 


The fifth moment is 


(60) 


, L-l / M-l k 

s 5 -Z Z { Z n x 

' i=0 V j=0 7=1 j,T 


(61) 


or 


, , L-l 

S 5 =i Z 


1=0 


1 + 5 (M-l) + 30 (M-l) (M-2) 


M-l k 

Z n x 

j=o y=l 


i+j.r 


M-3 M-2 M-l k 

+ (lOM)3l Z Z Z n x. . x. x. 


j=o v=j+l £=v+l 7=1 


i+j.Y i+v,7 i+£,y 


M-5 M-4 M-3. M-2 M-l k 

+ 51 Z Z Z Z Z 11 x i+i 


j=o v=j+l £=v+l f=£+l g=f+l 7=1 


X X X X 

i+j,y i+v,7 i+f,7 i+f,7 i+g,7 


(62) 
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Now distributing the sum over i and simplifying. 


L-l M-l k m-i k y 

Z Z n x i+i y = Z n Z Vi'-y = M( ' l) (63) 

i=o j=o 7=1 J j=o 7=1 i=o 1 


M-l k 


L -1 


L-l M-3 M -2 M-l k 

Z Z Z Z n x i+ j 

i=o j=o v=j+l f=v+l 7=1 J 


i+j»r X i+v,7 X i+S,7 


L -1 


M-3 M-2 M-l k V 

= 2 1 S 1 i Vi. Vr ^.r 

j=o v=j+l £=v+l 7=1 i=o J,T ,T 6,7 


= Z c n L (- 1 ) 1 

co , y 

co=l y=l J 


L-l M-5 M-4 M-3 M-2 M-l k 

i ?0 h vjl J=v + 1 «+l Jf+1 A X ‘ + j>'>' Xl+V -^ X ‘ +f ^ X ‘ +£ -->' X ‘ + 8.V 


M-5 M-4 M-3 M-2 M-l \ k 

v=^+ 1 fUl f| + l O-W A X ‘ + )^ Xl+v *->' X * +f -'>' ’W 


2 w . 

Z G n L (-1 ) 1 

, <o , y 

co=l L y=l 
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where M < Lj for all y, and G is the number of combinations of j* v, £, i, 

*• co 

^possible for which the set w (with w members) of the k 
component sequences satisfies 

( 66 ) 


Assuming that no common trinomials or pentanomials are satisfied for the 
characteristic equations of the maximum-length sequences comprising the 
sum sequence, equation (61) becomes 


X = X X X 

i+j,y i+v,y i+£,y i+g,y 



S 5 = [1 + 5 (M-l) + 30 (M-l) (M-2)] M(-l) k 

L 


10M31 


- X B (-i) k X B L (- i ) k-1 

\ 3 / r + yy 


T [(") - F J ( - 1)k + | F r L v ( - 1)k ' 1 ' <67) 


where M<L for all y, B is the number of trinomials of order M-l or less 

y th 7 

that contains the y -sequence characteristic polynomial as a factor, and F 

th 

is the number of pentanomials of order M-l or less that contains the y 
sequence polynomial as a factor. 

If M«L, 


X B L (- 1 

5 ^yy 

S « 10M3I 

c k 

n L 


+ 5! 


I W-D 

y=l r Y 


n l 
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is an approximation for the fifth-central moment. 


TEST ALGOR ITHM FOR PSEUDORANDOMNESS 
OF DIGITAL SEQUENCES 


The calculations performed previously for the first five central 
moments are given in Table 2. 


TABLE 2. APPROXIMATIONS FOR THE FIRST FIVE CENTRAL MOMENTS 



L = II L , and when no common trinomials, quadrinomials, and pentanomials 
y=l Y 

of order M-l or less contain the sequence characteristic polynomials as factors 
for each maximum-length sequence comprising the sum sequence. 
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A test involving the formation of a weighted sum of the difference in the 
first N moments for the sum sequence and a random sequence may be used to 
evaluate sequences from sum generators. This can be expressed as 


N / . . \ 

T(M) =• Yj W. IS 1 - s 1 
. , 1 V c cr 

i=l \ / 


(69) 


where S * is the i^-central moment for weights of M-tuples from a random 
cr ,, 

i th 

sequence, W. is a weighting factor, and S is the i -central moment for 
X c 

weights of M-tuples from a pseudorandom sequence. 

Using the results from Table 2, equation (69) becomes 


T(M) = 



W„ B 31 + W, 41 E 
3 y 4 y 


+ W '(10M31 B + 51 F ) 
5 y y' 


L (-1) 

y 


k-l 


(70) 


For a particular selection of the weighting functions, the smaller the value of 
T(M), the better the sequence approximates a random sequence. 

The weighting values can be selected to emphasize a particular aspect 
of the distribution of M -tuple weights. For example, the term 


w , l 


y=l l 


B L (-1) 

yy 


k-l 


n l 
Tl 1 


indicates the relative symmetry or skewing of the distribution. The term 
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(10M31 B + 5l F ) L (-l) k_1 

y y y 



indicates skewing of the distribution with more emphasis on the shape of the 
distribution of M-tuple weights beyond the variance of the distribution. The 
term 


k r 


W. 


y=l 


4! Ey (-1) 


k-1 



2 

indicates the kurtosis of the distribution of M-tuple weights. Assuming W > 0, 

3 4 

positive values of this term indicate a leptokurtic distribution, and negative 

values of the term indicate a platykurtic distribution. If k is odd the distribu- 
tion is leptokurtic, and if k is even the distribution is platykurtic. 

A computer algorithm for evaluating the sequence test parameter, 

T(M) , has been developed. The algorithm calculates B^, the number of tri- 
nomials of order M-l or less that contains the y^-sequence characteristic 
polynomial or a factor; E^, the number of quadrinomials of order M-l or less 

-that contains the y tk -sequence characteristic polynomial or a factor; and F_. 

th ^ 

the number of pentanomials of order M-l or less that contains the y -sequence 
characteristic polynomial or a factor. 


Lindholm developed an efficient algorithm for calculating B , and the 


algorithms developed for E and F 
method [7, p. 572]. 


are essentially extensions of Lindholm’s 


2. The relative flatness or peakedness. 

3. Flat compared with normal distribution. 
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If a sequence is generated by an n-stage register, and the sequence is 
a maximum-length type, then any 2n-l digits of the sequence can define the 
particular stages that contribute to the feedback. This is equivalent to solving 
n-1 simultaneous equations, since for a maximum -length sequence the last 
stage is always fed back. If the sequence characteristic polynomial is a 
factor of a trinomial of the form 


g(x) = 1 + x d - C + 


d 

X 


( 71 ) 


then the sequence satisfies the recursive relation 


x. x. 

1-C 1-1 


( 72 ) 


when the sequence is from the set {-1, +l} . 

One particular content vector in a maximum-length sequence is 

= -1, x = 1, x = 1, . . . x „ = 1, x = -1. Using this content vector as a 
0 1 2 n-1 n 

starting point, the next M-l content vectors are calculated using the sequence 

recursive relation. Then M + n digits of the sequence are known and can be 

represented as 


x 


O’ 


V 


X M + n - 1 


= - 1 , 1 , 1 , 1 , 


-l,x 


n+1’ n + 2* 


x. 


M + n - 1 


( 73 ) 


Because the tuple 


(x^» X^, • . . X Q ' j) — (l, 1 , . . . l) 


( 74 ) 
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and if the sequence characteristic polynomial is a factor of 

- , d-c d . . 

1 + x + x , (75) 

then the tuple 

(x j±1 , j i ■ « x. . ) (x ,, X , . . . X ) 

d + 1 d+2 d+n-l /v c + l c + 2 c + n - r 

= ( 1 , 1 , 1 , ... 1 , 1 ) . ( 76 ) 

If X is a vector representing the tuple (X, * x, x , ,) , and 

d v d+ld+2 c + n-1 7 

similarly for X , the relation 


* * 

X j X = I 

d c 


( 77 ) 


expresses the condition of equation ( 71) , where I is the identity matrix of order 

n-1, and X and X are (n-l) by (n-l) matrices with the elements of X , and 
d c d 

X , respectively, on the main diagonal and all other elements equal to zero, 
c 

Extending this procedure to quadrinomials and pentanomials that contain the 
sequence characteristic polynomial as factors, the required vector relations 
are 


* * * 
X J X X 
dee 


I 


and 


* * * 
X, X x 
d c 


e 
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( 78 ) 


( 79 ) 



By finding tuples for which these equations hold, using the first M-l 
digits after the content vector (l, 1, 1, ... -l) , all trinomials, quadri- 
nomials, and pentanomials that contain the sequence characteristic polynomial 
as a factor are yielded. 


The computer program POLTE 1 was written to solve for the vector 
relations in equations (76), (77), and (78). The results from this program 
can be used to evaluate T(M) from equation (69) . The procedure is as 
follows: 

1. Select M, the size of the M-tuple. 


2. Select k' sequences to form the sum sequence. 

3. Using the computerized algorithm, calculate B^, E^, and F^. 

4. Select the set of weightings, W., depending on the characteristics 
of the distribution of M-tuple weights that are critical. 

5. Evaluate T(M) from equation (69) . 

This procedure can be used to evaluate candidate designs for pseudo- 
random sequence generators of the type shown in Figure 1, 


Appendix A contains the results of the computer evaluation for a maxi- 
mum-length sequence with 11 memory elements and a hybrid-sum sequence 
with 11 memory elements. Appendix B is a flow diagram of the program 
POLTE 1. 


An indication of the reduction in the amount of calculations required to 
evaluate the statistics of a filtered hybrid-sum sequence as compared with a 
filtered maximum-length sequence can be determined as follows. The limit 
ratio of the number of pseudorandom sequences statistically evaluated com- 
pared with the amount of calculations required is 


R = 


k 

n 2 

i=l n 


n. 

l 


n. 


K 

l i l 

i=l n 


(80) 
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For k=l the ratio is unity, but as k increases the ratio tends to increase as 
previously illustrated. This means the hybrid-sum sequence generator 
configurations potentially can provide many pseudorandom digital sequences 
with a minimum number of calculations required. 

As an example of the potential increase in computational efficiency, 
using the hybrid-sum approach, assume that the computer algorithm was 
efficient enough so that each moment could be calculated in 1 second. It would 
then require more than 12 days of computer time to completely analyze the 
statistics of all possible maximum-length sequences from a 23 -stage register. 
If, however, the sequence group is established from the hybrid sum of 
sequences from 11- and 12-stage registers, the analysis of sequences, which 
are approximately 99.9 percent as long as the maximum- length sequences 
from the 23 -stage register, can be accomplished at a rate 120 times faster 
than the analysis in the maximum-length case. The more maximum-length 
sequences that form the hybrid-sum sequence, the greater the efficiency in 
forming the sequence in this manner. 

As an illustration of the theory presented, a comparison is made of the 
statistics of a filtered maximum-length sequence from the 11-stage generator 
and a filtered hybrid-sum sequence from a 5- and 6-stage generator. The 
filter impulse-response length is assumed to be 20 digital-clock periods. The 
11 -stage maximum-length sequence generator is described by the polynomial ,, 
(ll, 9, 0) and the hybrid-sum generator by the pair of polynomials (5, 2, 0) 
and (6, 1,0). 

An evaluation of equation (69) with , 

W equal to 1, 
equal to 0, 

and 

W_ equal to 0 
5 

gives an indication of the skewing of the amplitude distribution of filtered 
pseudorandom sequences. Using the described method of calculation, this 
parameter is evaluated as follows: 
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( 81 ) 


For filtered maximum-length sequence (ll, 9, 0), 

T(M=20, W = 1, W = 0, W = 0) = 54, 
o o 

indicating dominant positive skewing. 

For filtered hybrid-sum sequence (5, 2, 0) + (6, 1, 0), 

T(M=20, W = 1, W = 0, W = 0) = -8, (82) 

O 4 0 

indicating slight negative skewing. 

A computer program was written to evaluate the distribution of weights 
of the filtered sequence for both the maximum-length sequence and the hybrid- 
sum sequence directly. Figure 5 is the result for the maximum-length 
sequence. The distribution skews to the positive side and is a poor approxi- 
mation to the normal distribution. Figure 6 is the result for the hybrid-sum 
sequence and shows very little skewing tendency. 


APPLICATION TO COMMUNICATIONS SYSTEMS TESTING 

A diagram of an arrangement for communications system evaluation is 
presented in Figure 7. This test configuration includes a pseudorandom 
signal generator that supplies the test waveforms that simulate the waveforms 
many communications systems are required to process in normal applications. 
The evaluation scheme establishes as a system quality factor the root-mean- 
square (rms) value of the difference between system output and input. The 
test configuration also includes gain- and delay-compensation facilities, which 
are required because the system quality factor must be minimized with respect 
to these two parameters. During processing by a communications system, 
signals may be attenuated, delayed, distorted, and random noise may be 
added. Of these possibilities, attenuation and uniform delay are not con- 
sidered to be errors, since they represent no deterioration in the waveform’s 
ability to convey information. Delay and gain are therefore adjusted to match 
the average system delay and average system attenuation. The test configu- 
ration measures waveform distortion introduced into the communications 
process. 

The mathematical expression for the mean-square-error is [ 1 ] 
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OUTPUT AMPLITUDE 


Figure 5. Computer simulation results for 
maximum-length sequence (ll, 9, 0) . 

.(*) = t™* o? f [e.(t-r) -Ge Q (t) -Gn Q (t)] 2 dt , (83) 

-T 

where e.(t-r) is the delayed system input, e (t) is the system output, G is the 




20HHdWH?dH 
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OUTPUT AMPLITUDE 


Figure 6. Computer simulation results for 
hybrid-sum sequence (5, 2, 0) + (6, 1, 0) . 

gain-compensation factor, and n Q (t) is the noise added-in system. 

The value of the rms error measured as shown in Figure 7 depends 
upon the statistical and spectral properties of the systems input signal. There- 
fore, any specification relating to communications system performance as 
determined by the rms testing scheme must include a complete specification 
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AMPLIFIER 


»* 


Figure 7. System for communications network evaluation. 

of the test signal’s spectral density and amplitude probability distribution. 
Common requirements for system test signals are an approximately Gaussian 
amplitude probability density function and a signal spectral envelope with the 
shape of common filtering functions, such as the low -pass Butterworth or 
Chebychev functions of a particular order. 

The application of true random signals to simulate the input signals 
of the system is difficult because analog random signals cannot be delayed 
without some distortion. Also, the statistics of the noise signal must be 
stationary and the spectral characteristics must be well defined before the 
measurements can be assigned any real significance. Shepertycki's method 
of generating the system input signal and delayed input is shown in Figure 8. 
The advantage of this method is that the delay operation is performed digitally 
and virtually without distortion. The filters form the analog pseudorandom 
signal from the digital sequence. The transfer characteristics of the two 
filters must be closely matched to ensure proper operation. 

The hybrid-sum sequence can be the source sequence for the noise 
signal generation scheme of Figure 8. The primary reason for the hybrid-sum 
sequence approach is to form large groups of pseudorandom sequences with a 
minimum of computer evaluation required to establish an analysis of the 
statistics of the filtered sequences. It has been shown that groups of hybrid- 
sum sequences can be evaluated much more rapidly and efficiently than groups 
of maximum -length sequences. The hybrid-sum approach is especially useful 
if a multichannel communications system is to be evaluated by establishing the 
mean-square error parameter as given by equation (83) for each channel. 
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Figure 8. Method of generating test signals. 

A block diagram of the scheme for test signal generation and signal 
application for a multichannel communications system is shown in Figure 9. 
Eaph channel is evaluated on an rms basis, as illustrated in Figure 7. The 
test signals generated must be pseudoindependent in order for channel cross- 
talk to enter as the noise term, [n (t)I, in the error evaluation equation (83) . 



Figure 9. Test configuration for a multichannel communications system. 
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CONCLUSIONS 


The purpose of this study was to establish a procedure for designing 
and analyzing noise-like signals from filtered, pseudorandom digital sequences. 
An immediate application of these filtered sequences is as test signals for 
performance evaluation of a multichannel communications system. The method 
for multichannel communications system evaluation is essentially an extension 
of a scheme proposed by Shepertycki; he suggested using filtered maximum- 
length digital sequences as the source of the data signals. Previous studies 
performed by Gilson [3], White l6l, Lindholm [7l, andWainberg and Wolf 
[8] have shown that filtered maximum-length sequences do not necessarily 
have good statistical properties and therefore cannot be used indiscriminately 
in pseudorandom signal generators. 

This study proposes that test signals for this performance evaluation 
scheme be derived from the filtered modulo-two sum of several maximum- 
length digital sequences. Sequences generated in this manner are called 
hybrid-sum sequences. 

The theoretical analysis evaluates the statistics of filtered hybrid-sum 
sequences by calculating the higher moments of the weights of tuples from the 
sequence. This analysis is an extension of Lindholm's maximum-length se.- 
quences analysis to the hybrid-sum sequences. 

The theory developed by this study yields a procedure for evaluating 
large groups of potential code sequences with a level of efficiency not previ- 
ously possible; therefore, the results of the analysis may be applied not only 
to system testing but also to areas such as coding theory and modulation 
format design. 



APPENDIX A 

RESULTS OF COMPUTER PROGRAM POLTE 1 


This appendix contains example results of the computer program 
POLTE 1. The program was run for three irreducible polynomials of order 
11, 6, and 5. The results of POLTE 1 can be used to evaluate the statistics 
of filtered, pseudorandom digital sequences using equation (69) . 

The procedure for evaluating equation (69) using the results from 
POLTE 1 is as follows: For a given irreducible polynomial that generates a 
maximum-length sequence, the polynomial representation is printed in binary 
and octal form. For example, 


POLYNOMIAL 110000100000 

OCTAL 0103 

represents the polynomial 

6 5 , 
x + x 1. 

Following the polynomial octal form, the representations of the- tri- 
nomials, quadrinomials, and pentanomials that contain the characteristic 
polynomial as a factor are printed. For example, 

(7, 5, 1, 0) 

represents the polynomial 

7 5, 

X + X + X + 1. 

This quadrinomial contains the characteristic polynomial, 

x^ + x^ + x + 1 = (x^ + x^ + l) (x + l) (mod-2). (A-l) 

For a filter impulse-response period M, the parameters B^, E^, and 
F^ are respectively, the number of trinomials, quadrinomials, and penta- 
nomials of order M-l or less ' of the form 
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and 


P / d c , v 

x (x + x + 1) , 

p , d - c b x 

X (x + X + X + l) , 


p / d c b a , v 
x(x+x+x+x+l) 

til 

that contain the y -characteristic polynomial as a factor. In the above poly- 
nomials, p can range from 0 to M-l-d. Therefore, for each basic polynomial 
of the forms 

d c 
x + x + 1, 

deb 

x+x+x+1, 

or 

d c b a „ 
x + x + x + x + 1 

that contains the sequence characteristic polynomial as a factor, there are 
M-d product polynomials that also contain the sequence characteristic poly- 
nomial as a factor. 

The algorithm in the computer program POLTE 1 detects the number 
of basic polynomials that contains the sequence characteristic equation as a 
factor. The parameters B , E , and F are calculated from the results of 

y y y 

POLTE 1 by forming the sum 


B 

y 



Z ( M - d i) 

• -I ^ 


(A -2) 


E = 

y 


Z (M-d.)-f 

i= 1 


M-l 

Z 


i=l 


(M-i) 


(A-3) 
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and 


F y = Jj ( M ~ d i) > (A-4) 

where D , D 0 , and D are the number of trinomials, quadrinomials, and 
1 " 3 

pentanomials of order M-l or less that are detected for the y -component 

th 

sequence by the program POLTE 1. If the characteristic equation of the y 
sequence is a trinomial, f = 1; f = 0 otherwise. 

The program POLTE 1 was used to evaluate the parameter for 

maximum-length sequences (ll, 9, 0), (6, 1, 0) * and (5, 2, 0) , where M 
was equal to 20. The results are given in Table A-l. 

TABLE A -1. POLTE 1 FOR THREE 
MAXIMUM-LENGTH SEQUENCES 


SEQUENCE 

B 

Y 

(11, 9, 0) 

9 

(6, 1, 0) 

22 

(5, 2, 0) - 

38 


These values were previously used to evaluate the statistics of the 
maximum-length sequence (ll, 9, 0) and the hybrid-sum sequence (6, 1, 0) 
+ (5,2,0). 
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APPENDIX B 

FLOW DIAGRAM FOR THE PROGRAM POLTE 1 
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